Microscopic Theory and Simulation of Quantum- Well 
Intersubband Absorption 

Jianzhong Li and C. Z. Ning 

Center for Nanotechnology and NASA Advanced Supercomputing Division 
NASA Ames Research Center, Mail Stop N229-1, Moffett Field, CA 94035-1000 


A DCTD A 

XliJO JL JL 

We study the linear inter subband absorption spectra of a 15 nm In As quantum well using the intersubband 
semiconductor Bloch equations with a three-subband model and a constant dephasing rate. We demonstrate 
the evolution of intersubband absorption spectral line shape as a function of temperature and electron density. 
Through a detailed examination of various contributions, such as the phase space filling effects, the Coulomb 
many- body effects and the nonparabolicity effect, we illuminate the underlying physics that shapes the spectra. 
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1. INTRODUCTION 

Study of semiconductor optoelectronic devices in the infrared regime attracts great attention owing to their 
applications in communications, sensing, imaging, and “finger printing” molecules. Current efforts of extending 
the spectral range of available light sources and detectors toward long micrometer wavelength regime involve 
very diverse device concepts and designs, such as quantum cascade lasers and type II interband lasers. 1 One of 
the promising device designs for achieving longwave infrared, or far-infrared lasing is based on optically pumped 
intersubband transitions (ISBTs). 2 The advantages of this concept include large dipole matrix elements and 
ultrafast dynamic response. Such design could be made even more attractive if material systems with deeper 
quantum wells (QWs) are adopted, such as the InGaAs/AlAsSb QWs, since then the pumping source for such 
systems could be another compact semiconductor diode laser. 3 In this respect, the antimonide-based QWs 
offer an interesting alternative, as the conduction band offset of InAs/AlSb heterostructure is as deep as 2 eV, 
providing more design flexibilities. 

We have investigated effects of bandstructure, Coulomb interaction, electron-longitudinal optical (LO) phonon 
interaction, and other material-related issues on ISBTs within a two-subband model. 4 ” 7 The analysis is based 
on density matrix formalism. In this paper, linear intersubband absorption spectra in a 15 nm InAs quantum 
well are studied in the same framework — the intersubband semiconductor Bloch equations approach — but under 
a three- subband model and dephasing rate approximation. We demonstrate the evolution of intersubband 
absorption fine shape as temperature and carrier density are changed. Through a detailed examination of 
various contributions, such as phase space filling effects and Coulomb interaction-induced many-body effects, we 
elucidate the underlying physics that shapes the spectra. 

Theoretical treatment of intersubband transitions typically involves two steps. The first step is to compute 
the bandstructure, while the second step is to formulate light-semiconductor interaction on the basis of this 
obtained bandstructure. The first step is sometimes called determination of the ground state, while the second 
step treats the excitations of the system by external light field. In the present model study, we model the band- 
structure with different effective masses for the three parabolic subbands, whereas the density matrix theory 8,9 
is adopted to treat the light-semiconductor heterostructure interaction. By explicitly formulating the two-point 
correlation functions, such as the intersubband polarization, the set of the so-called intersubband semiconductor 
Bloch equations (iSBEs) is derived. Within this approach, it is known that the exchange interaction leads to 
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a self-energy renormalization (exchange self-energy, or XSE) to the single particle energy and a nonlocal ver- 
tex (excitonlike in the case of interband transitions) term that couples the other intersubband polarizations to 
each individual fc -transition, whereas the direct (Hartree) interaction contributes the so-called depolarization 
field term that describes the dynamic screening of the charge carriers to the external exciting field, which is of 
collective nature and modifies the local field that each individual carrier feels. This paper detailedly examines 
the case of a 15 nm In As quantum well within a three- subband model and illustrates how different contributions 
help shape the ISBTs as a fraction of temperature and carrier density. 

This paper is organized as follows: In the second section, we summarize our theoretical considerations with 
a minimum set of equations; we present simulated ISBT spectra for the In As quantum well in Sec. 3; and then 
we conclude the paper with a summary. 


2. THEORETICAL CONSIDERATIONS 

As mentioned above the microscopic theory consists of two steps: the static description of the subband disper- 
sions and kinetic description for the intersubband optical transitions. We describe the energy dispersions in an 
approximate manner in this work, or more specifically, different effective masses are used for the three subbands. 
The focus of our work is on the kinetic description of the ISBTs. 

To treat the fight-semiconductor heterostructure interaction, we consider the following expectation values, 
of the bilinear combination of creation (c^ fc ) and annihilation (cnfc) operators among the same k (the in-plane 
electron wavevector) states within a three-subband model (subband labeled by index n = 1, 2, 3): ground subband 
population f\k = (c| fc cifc), lower excited subband population fok = (4fc C2fc )’ u PP er excited subband population 
fzk = (4fc C3fc )’ intersubband polarization p 12 (k) ~ ( c\ k C 2 k ) between subband 1 and 2, p 23 (k) = {c\ k czk) between 
subband 2 and 3, p u {k) = { c\ k czk ) between subband 1 and 3, and their corresponding Hermitian conjugates. 
Following the quantum kinetic approach, the semiclassical kinetic equations for the above dynamic variables axe 
derived as a limiting case. 9 The derivation is extensive and will not be given here. The resultant equations 
(only the ones for inter subband polarizations here), under linearization with respect to the incident fight field 
amplitude, some plausible assumptions, and the rotating wave approximation, are found as follows: 

[h {U + i 7p mn ) - - £ mk)} PT = (dr • Eo + eD (fmk - fnk) ~ £ f jk+g {V^ pt ~ > (1) 

where lj is the angular frequency of the incident light of amplitude Eo, 7™ n is the dephasing rate, <i™ n is the 
dipole matrix element, and l ^ m, l ^ n. fmk is taken as the Fermi distribution function in the linear absorption 
calculation and p™ is the amplitude of the inter subband polarization p mn (k). The renormalized single particle 
energy ( e m k ) by the Coulomb self-energy and the local field correction term {e k n ) are, respectively, given by 


' • ( 2 ) 
l,q 

= - X K jml pk +q + X v o n3lm pH ■ ( 3 ) 

The local field correction consists of a Fock (first) term that gives rise to a type of collective excitation called the 
repellon and a Hartree (second) term that stimulates another type of collective excitation called the inter subband 
plasmons (ISPs). 6,10 We mention that the first term is responsible for the Fermi-edge singularity effect, whereas 
the remaining term leads to the depolarization effect, as described detailedly in our earlier work. Their effects 
have been known from previous studies and play the same roles as well within the present three- subband model. 
The Coulomb matrix elements (V^ mh s) are defined as in Ref.. 5 The static single plasmon-pole approximation 8 
has been used for screening the exchange interaction by the intrasubband processes in the present simulation 
work. The iSBEs are solved numerically by a matrix inversion for the intersubband polarization functions. We 
note that the last term in Eq. (1) is unimportant in this work. 

The linear absorption coefficient is defined by 
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with n{(jj) being the background index of refraction (slowly varying in frequency), e(w) the dielectric function, 
and c the speed of light in vacuo . The susceptibility, x(^)> ' 1S given by 

xH = P/eoEo , (5) 

P = 2S/[(27r) 2 V] J2 [ dk {C n »r • (6) 

m/n 

where P is the total intersubband polarization, sq is the electric constant, S is the QW area, and V is the QW 
volume. V = WS, and W is the QW thickness. Finally, the absorbance of the semiconductor heterostructure is 
given by 2 Wa(cu) (per bounce at full TM polarization). 

3. NUMERICAL RESULTS 

Table I lists the bandstructure-related parameter values that were used in the numerical simulation. All the 
values were obtained using the spurious-state-free 8-band k • p Hamiltonian under the envelope function approx- 
imation, 5,11 with the exception of the depolarization factors. Their values are the quantum box results of the 
same QW thickness. In order for data integrity, the calculated subband populations are also presented here in 
Table II, which are auxiliary to the understanding of the simulated spectra, as shown later. 


Table I. Parameters used in simulations 


subband effective mass 

Tfli (77i e ) 

0.0336 (1) 

0.049 (2) 

0.0715 (3) 

subband separation 

Eij (meV) 

131.065 (12) 

150.92 (23) 

281.985 (13) 

dipole matrix element 

dP (e-A) 

30.0 (12) 

35.0 (23) 

0.0 (13) 

depolarization factor 

D jlmn ( W ) 

0.1222 (1122) 

0.1146 (2233) 

0.0343 (1133) 0.1101 (1232) 


* m e : free electron mass; e: absolute electron charge. Numbers. in parenthesis after entry values indicate subband 
indices. See Ref. 10 for definition of the depolarization factor. 


The effective masses were found after a least square fitting to the k p subbands; the subband separation is the 
energy difference at the T point (k — 0) between the k p subbands; the dipole matrix elements are found to be 
weakly dependent on the wavevector so that their values at the T point were used in the simulation. We mention 
that the depolarization factor is a measure of the contributing strength of the Hartree term to the local field , 
which in turn determines the strength of the ISPs. 6,10 Note that a value of the order of 0.1, in the unit of the 
QW thickness (W), means a rather strong depolarization effect; the total strength of the depolarization effect is 
proportional to the product of the QW thickness and the density difference between the subbands. Furthermore, 
coherent Coulomb effects due to the coupling of ISPs associated with individual ISBTs are important, and the 
results will be presented elsewhere. 


Table II. Calculated subband populations 


electron density (10 12 cm ^) 

temperature (K) 

m 

n 2 

nz 

1.00 

12 

1.00 

0.00 

0.00 


80 

1.00 

2.41 xlO- 5 

1.09 XlO- 14 


300 

9.58 xlO" 1 

4.16X10- 2 

1.84 xlO" 4 

2.50 

12 

2.11 

3.92 x nr 1 

0.00 


80 

2.10 

3.95 xlO -1 

9.91x10- 10 


300 

2.00 

4.96xl0 _1 

3.49xl0 -3 

5.00 

12 

3.13 

1.87 

0.00 


80 

3.13 

1.87 

3.77xl0- 5 


300 

3.09 

1.84 

6.78 xl0“ 2 

7.50 

12 

4.06 

3.23 

2-lOxlO -1 


80 

4.04 

3.21 

2.49X10- 1 


300 

3.93 

3.06 

5.10X10- 1 


rii: Electron density in Subband i. 




Shown in Fig. 1 are the energy dispersions of the considered three conduction subbands (left panel) and the 
single-particle absorption spectrum (right panel) of the 15 nm In As QW. Owing to a nonparabohcity of the bulk 
InAs conduction band, dispersions for the three subbands are different; the subband separations are fc-dependent. 
The rather strong nonparabolicity of InAs leads to a large value range of the intersubband separation. In a single- 
particle picture, this introduces an inhomogeneous broadening and an accompanying low frequency tailing to 
the absorption spectrum. All these are reflected in the figure. Note that we did not consider the temperature 
dependence of the subband strucuture in this model study. Also shown in the figure, as a comparison, is the case 
of vanishing nonparabohcity (dotted curves with upward arrows): The spectrum consists of sharp resonances 
without any inhomogeneous broadening. 
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Figure 1 . Schematic of subband dispersions (left panel) and single-particle intersubband absorption spectrum for a 15 nm 
InAs quantum well. The decreasing intersubband separation with increasing k represents the result of nonparabolicity in 
InAs conduction band, which is consequently and unambiguously reflected by the broadened absorption spectrum. 


As out fined in Sec. 2, the matrix inversion method was used to solve Eq. (1) for the intersubband polarization 
functions. Then the absorbance was calculated using these functions according to Eqs. (4)-(6) and the results 
are presented in Figs. 2-4. We first show how the absorbance changes as a function of temperature at different 
electron densities in Fig. 2. Since the two-subband results have been understood rather well, 6, 10 we thus 
emphasize density range where a three-subband model is warranted, that is when the lower excited subband 2 
starts to be populated. We have chosen four densities, as denoted in the figures as well as in Table II. At the 
lowest one, the second subband is only thermally populated at 300 K, as witnessed by the appearance of the weak 
resonance from the lower excited subband 2 to the upper excited subband 3 (2— >-3) near 150 meV. Otherwise, 
only a single, inhomogeneous ly broadened resonance (1— >2) is present. Expectedly, it is further weakened and 


broadened thermally, as electrons increasingly populate higher energy states that have smaller transition energies, 
as displayed in Fig. 1. Increase in electron density strengthens resonance 2—^3 because the amplitude of the 
resonance is roughly proportional to the density difference of the two subbands, n 2 — n$. This dependence on 
the density difference also provides an explanation for the reduction in resonance strength as temperature rises 
at the four densities, which tends to deminish the difference. This is a phase space filling effect. Its another 
manifestation is seen in the red shift of resonance 1— >2 — more distinctly in Fig. 3 — as a function of density. 
On the other hand, note that at certain density, 5xl0 12 cm -2 here, resonance 2— >3 is temperature insensitive, 
which is a Coulomb collective effect. As we understand it now, ISBT is a collective phenomenon associated with 
the inter subband plasmon under normal circumstances, which could be a robust response, depending upon how 
strong the plasmon couples to the decay channels, such as through Landau damping. 12 
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Figure 2. Absorbance evolution as a function of temperature at different electron densities. At low density, only the 
ground subband 1 is occupied and thus only single resonance (1— >2) is observed at low temperature. Thermal population 
of the lowere excited subband 2 (see Table II) is revealed by the presence of a second resonance (2— >3). This resonance 
grows as density increases. Accordingly, the first resonance is weakened and redshifted, as a consequence of phase space 
filling. Note that at certain density, like 5xl0 12 cm -2 , the second resonance is temperature insensitive. 

The same results as in Fig. 2 are further presented as a function of electron density at different temperatures 
in Fig. 3 to illustrate all the effects from a different perspective. The phase space filling effects are more directly 
demonstrated in this case. First, at this particular QW thickness (15 nm), the first resonance line shape at 
the lowest chosen density of lxlO 12 cm -2 is single-particle-like, or a spectrum reflecting the constant 2D joint 
density of states as a result of the interplay of collective excitations. 6 Increase in electron density has two 
consequences: starting populating the higher subband(s) and thus redistributing the oscillator strengths. As 
seen from resonance l-»2, a larger density broadens its line shape further, tilts it toward the high-frequency side, 
and redshifts the spectrum. All these are indications that the ISPs play a more dominant role in the interplay, 
which is expected for higher density cases. Then the resonance is gradually diminished as further increase in 
electron density actually decreases the density difference of the two subbands, n\ — n 2 . Also, more redshift 








of the spectrum is observed. On the other hand, resonance 2— >3 starts to appear below a density of 2xl0 12 
cm" 2 and takes a symmetric line shape at low temperature. This is changed as either the temperature or the 
density increases, as seen in the figure. Within the present three-subband model, further increase in electron 
density first blueshifts the spectrum before the upper subband begins to be populated, because the self-energy 
renormalization of the populated subband increases the effective intersub band separation — similar to the two- 
subband model results. 6 However, as expected, continuing increasing the electron density will start to populate 
the upper subband 3 (see Table II). Then, the oscillator strength will eventually decrease and the spectrum will 
redshift, just like resonance 1— »2. Ultimately, we need to consider resonance(s) from subband 3 to even higher 
one(s). Note that we have not explicitly presented simulation results for resonance 1 — ^3, which lies at higher 
energy that is beyond the scope of our interested spectral range, but its physics is fundamentally the same as 
what we have discussed with regard to the other two resonances. 
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Figure 3. Same absorbance data as in Fig. 2, but presented as a function of electron density at different temperatures. 
Clearly seen is the phase space filling effects: Occupancy of the lower excited subband 2 introducing the second resonance 
(2— >-3) and adding to its oscillator strength as density increases; the enhancement of the plasmon peak at low temperature 
(12 K); redshifting both resonances as lower k states are occupied. Note the anomalous decrease at 300 K of the second 
resonance — a reflection of the robustness of the plasmon excitation: A phenomenon also observed by Warburton et aL. 12 

What has not been touched upon till now is the nonparabolicity effect, and this is the focus of Fig. 4. To 
demonstrate that, we choose an extreme case, that is to set the effective mass (m3) of the upper excited subband 
equal to that (m2) of the lower excited subband. In other words, there is no nonparabolicity effect, to the 
lowest order, for resonance 2— >3 in the present case. Furthermore, the intersubband separation E 2 3 is set to 
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173 meV, which bears no particular meaning other than to avoid the subject of the coupling of intersubband 
plasmons, which is the topic of a future paper. Not surprisingly, resonance 1 — >2 behaves as what have been shown 
previously. Furthermore, as expected, we observe narrow and symmetric Lorentzian line shape for resonance 
2-*3. Two features are worth noting: (i) an anomalous temperature dependence of the resonance at 2xl0 12 
cm~ 2 and (ii) its temperature insensitivity at higher density. On the first feature, the resonance is strongly 
enhanced when temperature rises. This turns out to be a manifestation of the phase space filling effects: At 
low temperature, the second subband is populated with a rather small electron density. As the temperature 
increases, more electrons are thermally excited from the ground subband into the second one. As a result of this, 
the number of electrons available for absorbing incident photons at resonance 2— >3 increases as well. Thanks to 
zero nonparabolicity, the temperature rise does not broaden the fine shape, in contrast to a large nonparabolicity 
case as discussed already. Therefore, the anomalous temperature dependence of the resonance at the density of 
2xl0 12 cm -2 is observed. Following the same line of thinking, the second feature of temperature insensitivity 
at higher density is easily understood. This also agrees with Warburton et al. 12 Finally, the blueshift of the 
resonance with increase in electron density reflects an effect of the Hartree contribution in Eq. (3). It is often 
called the depolarization shift in literature 13, 14 as, in the present case of zero nonparabolicity, it acts similar to 
a depolarization field. 
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Figure 4. Absorbance evolution similar to Fig. 2, but depicting the nonparabolicity effect as the effective mass of the 
upper excited subband was set to equal that of the lower excited subband. The intersubband separation E23 is set to 
173 meV. Two features are worth noting regarding the second resonance: (i) anomalous temperature dependence at 
2xl0 12 cm' 2 and (ii) its temperature insensitivity at higher density. 




4. SUMMARY 

In conclusion, we present a microscopic theoretical approach to intersubband optical resonances, similar to the 
semiconductor Bloch equations approach to interband transitions. This approach is applied to a 15 nm In As 
quantum well case that is representative of materials with large conduction band nonparabolicity. The evolution 
of intersubband resonances as a function of the temperature and the electron density is presented and discussed 
in detail. We demonstrate the importance of the phase space filling effects and Coulomb many-body effects in 
understanding the underlying physics of inter subband resonances. Finally, we show that the nonparabolicity 
effect is instrumental in shaping the spectra of intersubband resonances. 
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ABSTR A 

We study the linear intersubband absorption spectra of a 15 rnn InAs quantum well using the intersubband 
semiconductor Bloch equations with a three-subband model and a constant dephasing rate. We demonstrate 
the evolution of intersubband absorption spectral line shape as a function of temperature and electron density. 
Through a detailed examination of various contributions, such as the phase space filling effects, the Coulomb 
many-body effects and the nonparabolicitv effect, we illuminate the underlying physics that shapes the spectra. 
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1. INTRODUCTION 

Study of semiconductor optoelectronic devices in the infrared regime attracts great attention owing to their 
applications in communications, sensing, imaging, and “finger printing” molecules. Current efforts of extending 
the spectral range of available light sources and detectors toward long micrometer wavelength regime involve 
very diverse device concepts and designs, such as quantum cascade lasers and type II interband lasers. 1 One of 
the promising device designs for achieving longwave infrared, or far-infrared lasing is based on optically pumped 
intersubband transitions (ISBTs). 2 The advantages of this concept include large dipole matrix elements and 
ultrafast dynamic response. Such design could be made even more attractive if material systems with deeper 
quantum wells (QWs) are adopted, such as the InGaAs/AlAsSb QWs, since then the pumping source for such 
systems could be another compact semiconductor diode laser. 3 In this respect, the antimonide-based QWs 
offer ail interesting alternative, as the conduction band offset of InAs/AlSb heterostructure is as deep as 2 eV, 
providing more design flexibilities. 

We have investigated effects of bandstructure, Coulomb interaction, electron-longitudinal optical (LO) phonon 
interaction, and other material-related issues on ISBTs within a two-subband model. 4 " 7 The analysis is based 
on density matrix formalism. In this paper, linear intersubband absorption spectra in a 15 nm InAs quantum 
well are studied in the same framework — the inter subband semiconductor Bloch equations approach — but under 
a three-subband model and dephasing rate approximation. We demonstrate the evolution of intersubband 
absorption hue shape as temperature and carrier density are changed. Through a detailed examination of 
various contributions, such as phase space filling effects and Coulomb interaction-induced many-body effects, we 
elucidate the underlying physics that shapes the spectra. 

Theoretical treatment of intersubband transitions typically involves two steps. The first step is to compute 
the bandstructure, while the second step is to formulate light- semi conductor interaction on the basis of this 
obtained bandstructure. The first step is sometimes called determination of the ground state, while the second 
step treats the excitations of the system by external light field. In the present model study, we model the band- 
structure with different effective masses for the three parabolic subbands, whereas the density matrix theory 8,9 
is adopted to treat the light-semiconductor heterostructure interaction. By explicitly formulating the two-point 
correlation functions, such as the intersubband polarization, the set of the so-called intersubband semiconductor 
Bloch equations (iSBEs) is derived. Within this approach, it is known that the exchange interaction leads to 
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a self-energy renormalization (exchange self-energy, or XSE) to the single particle energy and a nonlocal ver- 
tex (excitonlike in the case of interband transitions) term that couples the other intersubband polarizations to 
each individual fc-tr ansi t ion, whereas the direct (Hartree) interaction contributes the so-called depolarization 
field term that describes the dynamic screening of the charge carriers to the external exciting field, which is of 
collective nature find modifies the local field that each individual carrier feels. This paper detailedly examines 
the case of a 15 nm InAs quantum well within a three- subband model and illustrates how different contributions 
help shape the ISBTs as a fraction of temperature and carrier density. 

This paper is organized as follows: In the second section, we summarize our theoretical considerations with 
a minimum set of equations; we present simulated ISBT spectra for the InAs quantum well in Sec. 3; and then 
we conclude the paper with a summary. 


2. THEORETICAL CONSIDERATIONS 

As mentioned above the microscopic theory consists of two steps: the static description of the subband disper- 
sions and kinetic description for the intersubband optical transitions. We describe the energy dispersions in an 
approximate maimer in this work, or more specificalty, different effective masses are used for the three subbands. 
The focus of our work is on the kinetic description of the ISBTs. 

To treat the light-semi conductor heterostructure interaction, we consider the following expectation values 
of the bilinear combination of creation (c^ fc ) and annihilation (c n k) operators among the same k (the in-plane 
electron wavevector) states within a three-subband model (subband labeled by index n = 1, 2, 3): ground subband 
population fik = (cj fc cifc), lower excited subband population fok = { c 2 k c ^ k )' u PP er excited subband population 
/ 3fc = {cl k czk), intersubband polarization p 12 (k) = (cj fc C 2 fc) between subband 1 and 2, p 23 (k) = (c^csfc) between 
subband 2 and 3, p 13 (k) = (c^ fc C 3 fe) between subband 1 and 3, and their corresponding Hermitian conjugates. 
Following the quantum kinetic approach, the semiclassical kinetic equations for the above dynamic variables are 
derived as a limiting case. 9 The derivation is extensive and wall not be given here. The resultant equations 
(only the ones for intersubband polarizations here), under linearization with respect to the incident light field 
amplitude, some plausible assumptions, and the rotating wave approximation, are found as follows: 

[n (w + *£“) - (s nk - e mk )\ P r = (dr • Eo + eJT) (fmk - U) - £ /,-*+, (V^pf - V™p£) , ( 1 ) 


w T here uj is the angular frequency of the incident light of amplitude Eo, 7™’ n is the dephasing rate, d™ n is the 
dipole matrix element, and l ^ m, I ^ n. fmk is taken as the Fermi distribution function in the linear absorption 
calculation and p™ a is the amplitude of the intersubband polarization p mn (k). The renormalized single particle 
energy (c m *.) by the Coulomb self-energy and the local field correction term (e™ n ) are, respectively, given by 

^ = i$L-'Z v ? bul f*+*’ ■ ( 2 ) 

=T = - X ^Pfc+«, + X V o 3lm ^‘ ■ ( 3 ) 

jrl^q j^l-.q 

The local field correction consists of a Fock (first) term that gives rise to a type of collective excitation called the 
repellon and a Hartree (second) term that stimulates another type of collective excitation called the intersubband 
plasmons (ISPs). 6,10 We mention that the first term is responsible for the Fermi-edge singularity effect, whereas 
the remaining term leads to the depolarization effect, as described detailedly in our earlier work. Their effects 
have been known from previous studies and play the same roles as well within the present three-subband model. 
The Coulomb matrix elements (V^ mh s) are defined as in Ref.. 5 The static single plasmon-pole approximation 8 
has been used for screening the exchange interaction by the intrasubband processes in the present simulation 
work. The iSBEs are solved numerically by a matrix inversion for the intersubband polarization functions. We 
note that the last term in Eq. (1) is unimportant in this work. 

The linear absorption coefficient is defined by 

= ;~7 , Em{e(u})} « — ^Im{x(u)} , 
son[uj)c n(u))c 


(4) 



with n(u?) being the background index of refraction (slowly varying in frequency), eta) the dielectric function, 
and c the speed of light in vacuo. The susceptibility. 1S given by 

x(u) = P/e 0 E 0 , (5) 

P = 2S/[(2tt) 2 V] Y, f dk {C}‘fT • (6) 

m^n ' 

where P is the total intersubband polarization, so is the electric constant, S is the QW area, and V is the QW 
volume. V = W S. and W is the QW thickness. Finally, the absorbance of the semiconductor heterostructure is 
given by 2Wa(u>) (per bounce at full TM polarization). 

3. NUMERICAL RESULTS 

Table I lists the bandstructure-related parameter values that were used in the numerical simulation. All the 
values were obtained using the spurious- state-free 8-band k • p Hamiltonian under the envelope function approx- 
imation, 0 ’ 11 with the exception of the depolarization factors. Their values are the quantum box results of the 
same QW thickness. In order for data integrity, the calculated subband populations are also presented here in 
Table II, which are auxiliary to the understanding of the simulated spectra, as shown later. 


Table I. Parameters used in simulations 


subband effective mass 

rm (m e ) 

0.0336 (1) 

0.049 (2) 

6.0715 (3) 

subband separation 

Eij (meV) 

131.065 (12) 

150.92 (23) 

281.985 (13) 

dipole matrix element 

d ij (e-A) 

30.0 (12) 

35.0 (23) 

0.0 (13) 

depolarization factor 

D jlmn 

0.1222 (1122) 

0.1146 (2233) 

0.0343 (1133) 0.1101 (1232) 


• m e : free electron mass; e: absolute electron charge. Numbers in parenthesis after entry values indicate sub band 
indices. See Ref. 10 for definition of the depolarization factor. 


The effective masses were found after a least square fitting to the kp subbands; the subband separation is the 
energy difference at the F point ( k — 0) between the k • p subbands; the dipole matrix elements are found to be 
weakly dependent on the wavevector so that their values at the F point were used in the simulation. We mention 
that the depolarization factor is a measure of the contributing strength of the Hartree term to the local field , 
which in turn determines the strength of the ISPs. 6,10 Note that a value of the order of 0.1, in the unit of the 
QW thickness (W), means a rather strong depolarization effect; the total strength of the depolarization effect is 
proportional to the product of the QW thickness and the density difference between the subbands. Furthermore, 
coherent Coulomb effects due to the coupling of ISPs associated with individual ISBTs are important, and the 
results will be presented elsewhere. 


Table II. Calculated subband populations 


electron density (10 i2 cm 2 ) 

temperature (K) 

ni 

77,2 

n 3 

1.00 

12 

1.00 

0.00 

0.00 


80 

1.00 

2.41xl0 -5 

1.09xl0~ 14 


300 

9.58xl0 _1 

4. 16xl0- 2 

1.84xl0 -4 

2.50 

12 

2.11 

3.92 xlO" 1 

0.00 


80 

2.10 

3.95 xlO- 1 

9.91xlO~ 10 


300 

2.00 

4.96 xlO -1 

3.49 xlO -3 

5.00 

12 

3.13 

1.87 

0.00 


80 

3.13 

1.87 

3.77xl0 -5 


800 

3.09 

1.84 

05 

Pj 

00 

X 

I— 1 

o 

1 

lO 

7.50 

12 

4.06 

3.23 

2.10XKT 1 


80 

4.04 

3.21 

2.49 xlO" 1 


300 

3.93 

3.06 

5-lOxlO- 1 


Tin Electron density in Subband i. 





Shown m Fig. 1 are the energy dispersions of the considered three conduction subbands (left panel) and the 
single-particle absorption spectrum (right panel) of the 15 mn InAs QW. Owing to a nonparabolicitv of the bulk 
InAs conduction band, dispersions for the three subbands are different: the subband separations are A;- dependent. 
The rather strong nonparabolicitv of InAs leads to a large value range of the intersubband separation. In a single- 
particle picture, this introduces an inhomogeneous broadening and an accompanying low frequency tailing to 
the absorption spectrum. All these are reflected in the figure. Note that we did not consider the temperature 
dependence of the subband strucuture in this model study. Also shown in the figure, as a comparison, is the case 
of vanishing nonparabolicitv (dotted curves with upward arrows): The spectrum consists of sharp resonances 
without any inhomogeneous broadening. 




k 


hv 


Figure 1 . Schematic of subband dispersions (left panel) and single- particle intersubband absorption spectrum for a 15 nm 
InAs quantum well. The decreasing intersubband separation with increasing k represents the result of nonparabolicity in 
InAs conduction band, which is consequently and unambiguously reflected by the broadened absorption spectrum. 

As outlined in Sec. 2, the matrix inversion method was used to solve Eq. (1) for the intersubband polarization 
functions. Then the absorbance was calculated using these functions according to Eqs. (4)-(6) and the results 
are presented in Figs. 2-4. We first show how the absorbance changes as a function of temperature at different 
electron densities in Fig. 2. Since the two-subband results have been understood rather well, 6, 10 we thus 
emphasize density range where a three-subband model is warranted, that is when the lower excited subband 2 
starts to be populated. We have chosen four densities, as denoted in the figures as well as in Table II. At the 
lowest one, the second subband is only thermally populated at 300 K, as witnessed by the appearance of the weak 
resonance from the lower excited subband 2 to the upper excited subband 3 (2—^3) near 150 meV. Otherwise, 
only a single, inhomogeneously broadened resonance (1— >2) is present. Expectedly, it is further weakened and 


broadened thermally, as electrons increasingly populate higher energy states that have smaller transition energies, 
as displayed in Fig. 1. Increase in electron density strengthens resonance 2—^3 because the amplitude of the 
resonance is roughly proportional to the density difference of the two subbands, no — This dependence on 
the density difference also provides an explanation for the reduction in resonance strength as temperature rises 
at the four densities, which tends to deminish the difference. This is a phase space filling effect. Its another 
manifestation is seen in the red shift of resonance 1-+2 — -more distinctly in Fig. 3 — as a function of density. 
On the other hand, note that at certain density, 5xl0 12 cm -2 here, resonance 2—^3 is temperature insensitive, 
which is a Coulomb collective effect- As we understand it now, ISBT is a collective phenomenon associated with 
the intersubband plasmon under normal circumstances, which could be a robust response, depending upon how 
strong the plasmon couples to the decay channels, such as through Landau damping. 12 
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Figure 2. Absorbance evolution as a function of temperature at different electron densities. At low denshy, only the 
ground subband 1 is occupied and thus only single resonance (1—^2) is observed at low temperature. Thermal population 
of the lowere excited subband 2 (see Table II) is revealed by the presence of a second resonance (2— »3). This resonance 
grows as densit}' increases. Accordingly, the first resonance is weakened and redshifted, as a consequence of phase space 
filling. Note that at certain density, like 5xlO L2 cm" 2 , the second resonance is temperature insensitive. 

The same results as in Fig. 2 are further presented as a function of electron density at different temperatures 
in Fig. 3 to illustrate all the effects from a different perspective. The phase space filling effects are more directly 
demonstrated in this case. First, at this particular QW thickness (15 nm), the first resonance line shape at 
the lowest chosen density of lxlQ 12 cm~ 2 is single-particle-like, or a spectrum reflecting the constant 2D joint 
density of states as a result of the interplay of collective excitations. 6 Increase in electron densit}' has two 
consequences: starting populating the higher subband(s) and thus redistributing the oscillator strengths. As 
seen from resonance l—>2, a larger density broadens its line shape further, tilts it toward the high-frequency side, 
and redshifts the spectrum. All these are indications that the ISPs play a more dominant role in the interplay, 
which is expected for higher density cases. Then the resonance is gradually diminished as further increase in 
electron density actually decreases the density difference of the two subbands, n\ — no* Also, more redshift 







of the spectrum is observed. On the other hand, resonance 2—5-3 starts to appear below a density of 2xl0 12 
cm" 2 and takes a symmetric line shape at low temperature. This is changed as either the temperature or the 
density increases, as seen in the figure. Within the present three-subband model, further increase in electron 
density first blueshifts the spectrum before the upper subband begins to be populated, because the self-energy 
renormalization of the populated subband increases the effective intersubband separation — similar to the two- 
subband model results. 6 However, as expected, continuing increasing the electron density will start to populate 
the upper subband 3 (see Table II). Then, the oscillator strength will eventually decrease and the spectrum will 
redshift, just like resonance 1— >2. Ultimately, we need to consider resonance(s) from snbband 3 to even higher 
one(s). Note that we have not explicitly presented simulation results for resonance 1^3, which lies at higher 
energy that is beyond the scope of our interested spectral range, but its physics is fundamentally the same as 
what we have discussed with regard to the other two resonances. 
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Figure 3. Same absorbance data as in Fig. 2, but presented as a function of electron density at different temperatures. 
Clearly seen is the phase space filling effects: Occupancy of the lower excited subband 2 introducing the second resonance 
(2— +3) and adding to its oscillator strength as density increases; the enhancement of the plasmon peak at low temperature 
(12 K); redshifting both resonances as lower k states are occupied. Note the anomalous decrease at 300 K of the second 
resonance — a reflection of the robustness of the plasmon excitation: A phenomenon also observed by Warburton et aL 12 

What has not been touched upon till now is the nonparabolicity effect, and this is the focus of Fig. 4. To 
demonstrate that, we choose an extreme case, that is to set the effective mass (m 3 ) of the upper excited subband 
equal to that (m 2 ) of the lower excited subband. In other words, there is no nonparabolicity effect, to the 
lowest order, for resonance 2^3 in the present case. Furthermore, the intersubband separation £33 is set to 




173 meV. which bears no particular meaning other than to avoid the subject of the coupling of intersubband 
plasmons, which is the topic of a future paper. Not surprisingly, resonance 1-+2 behaves as what have been shown 
previously. Furthermore, as expected, we observe narrow and symmetric Lorentzian line shape for resonance 
2-+3. Two features are worth noting: (i) an anomalous temperature dependence of the resonance at 2xl0 12 
cin -2 and (ii) its temperature insensitivity at higher density. On the first feature, the resonance is strongly 
enhanced when temperature rises. This turns out to be a manifestation of the phase space filling effects: At 
low temperature, the second subband is populated with a rather small electron density. As the temperature 
increases, more electrons are thermally excited from the ground subband into the second one. As a result of this, 
the number of electrons available for absorbing incident photons at resonance 2-^3 increases as well. Thanks to 
zero nonpar aboli city, the temperature rise does not broaden the line shape, in contrast to a large nonparabolicity 
case as discussed already. Therefore, the anomalous temperature dependence of the resonance at the density of 
2xl0 12 cm 2 is observed. Following the same line of thinking, the second feature of temperature insensitivity 
at higher density is easily understood. This also agrees with Warburton et al. 12 Finally, the blueshift of the 
resonance with increase in electron density reflects an effect of the Hartree contribution in Eq. (3). It is often 
called the depolarization shift in literature 13,14 as, in the present case of zero nonparabolicity, it acts similar to 
a depolarization field. 



Figure 4. Absorbance evolution similar to Fig. 2, but depicting the nonparabolicity effect as the effective mass of the 
upper excited subband was set to equal that of the lower excited subband. The intersubband separation E 23 is set to 
173 meV. Two features are worth noting regarding the second resonance: (i) anomalous temperature dependence at 
2xl0 12 cm -2 and (ii) its temperature insensitivity at higher density. 



4. SUMMARY 

In conclusion, we present a microscopic theoretical approach to intersubband optical resonances, similar to the 
semiconductor Bloch equations approach to interband transitions. This approach is applied to a 15 nm InAs 
quantum well case that is representative of materials with large conduction band nonpar abolicity. The evolution 
of intersubband resonances as a function of the temperature and the electron density is presented and discussed 
in detail. We demonstrate the importance of the phase space filling effects and Coulomb many-body effects in 
understanding the underlying physics of intersubband resonances. Finally, we show that the nonpar abolicity 
effect is instrumental in shaping the spectra of intersubband resonances. 
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